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Abstract 

The influence of a magnetic field on the tunnehng of an electron out of a confining 
plane is studied by a path integral method. We map this 3-d problem on to a 1-d one, 
and find that the tunneling is strongly affected by the field. Without a perpendicular 
field the tunneling at zero temperature can be completely suppressed by a large parallel 
field, but in the small parallel field and low temperature limit the tunneling rate is finite. 
An explicit formula is obtained in this case. A quantitative explanation without fitting 
parameter to a recent experiment is provided. 

PACS#s: 73.40.Gk; 73.20.Dx; 
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The tunneling of charged particles in the presence of a magnetic field has many peculiar 
properties. It occurs in the study of the escape of trapped electrons from a helium-vacuum 
interface [0 ^, the focus of the present paper. It is, however, a general phenomenon 
and is widely encountered in mesoscopic systems 0], in semiconductor heterostructures[|], ||, 
in disorder media[^, and in Quantum Hall Systems such as the tunneling between double 
layers|§] and the quasiparticle tunneling^]. The helium- vacuum-interface system is relatively 
simple, but it shares many features of the systems in Refs.[^ |^, ^ |^, |, ^, such as the 
influences of magnitude and direction of a magnetic field and of many-body correlations on 
the tunneling. This system provides us a well controlled test ground for our understanding 
of tunneling in higher dimensions. There is a considerable experimental study of the escape 
rate on it[]T], |^. In the case of no magnetic field the tunneling process in this system is 
adequately understood, which is effectively a 1-d problem even after the consideration of 
many-body correlations 0. However, the presence of a magnetic field requires a full 
3-d study of the tunneling process. It gives rise to new effects about the magnetic field 
and temperature dependence of the tunneling rate, as shown in a recent experiment which 
requires a clear and detailed theoretical explanation 0]. The purpose of the present paper 
is to give a formulation of the problem in this system and to present a detailed analytical 
study. A quantitative agreement without fitting parameter with the experiment of Ref . is 
obtained. Our results, as well as our path integral method which keeps the relevant degrees 
of freedom but eliminates irrelevant ones by integration, may have a wider applications in the 
systems of Refs.0, ||, |^, ^, |, |]. In fact, our Hamiltonian is identical to those, for example, 
in the study of magnetotransportation in semiconductor heterostructures[^]. 

Our main results are as follows. We consider a situation in 3-d that in one direction an 
electron has a metastable state, and there is no force in the remaining two directions. The 
tunneling out the metastable state is studied as a function of both magnitude and direction 
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of an externally applied constant magnetic field, and of temperature. We map this problem 
onto that of quantum dissipative tunneling. The in-plane motion acts as an environment 
to the tunneling in the perpendicular direction because of the parallel magnetic field. This 
environmental influence is found to be subohmic when the perpendicular component of the 
magnetic field Bj^ = 0, and superohmic when i?^ 7^ 0. From the obtained effective action we 
find that the tunneling is strongly influenced by the magnetic field. In particular, if B± = 
and the magnetic field B\\ parallel to the force free plane is large enough there is a complete 
suppression of tunneling rate at zero temperature, and the electron becomes dynamically 
localized. However, if B\\ is small, the tunneling rate out of the metastable state is finite. 
Perturbatively, the semiclassical action is proportional to Bp when 5^ = 0. In this case 
the calculation can be generalized to a finite but low temperature, where we find that the 
semiclassical action decreases linearly with temperature T. The last two results are in a 
remarkable agreement with the experiment of Ref . 0] . 

Now we present the calculation leading to the above results. For the convenience of 
calculation we shall first assume that the x- and y-directions have weak harmonic potentials. 
We take them to be zero, that is, ujq 0, in the end of the calculation. The Hamiltonian 
for an electron with mass m and charge e is then 

= i-[P - ^A(R)]^ + Imulix' + y') + V{z) , (1) 

with the vector potential A determined by V x A = B. Here the metastable point of 
the potential is taken at ^ = 0. For large enough z, V{z) = —Fz, and there is a barrier 
separating this region and z = 0. The external magnetic field is tilted: B = [B^O, Bj_). In 
accordance with the calculation of tunneling, the vector potential will be taken as 

A = (0,5^x-5||Z,0) , (2) 

and it can be shown that results are independent of the choice of gauge because of the 
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periodic boundary condition in the calculation of the tunneling rate. This is a 3-d tunneling 
problem because the magnetic field couples the motions in all directions. 
The tunneling is described by the Euclidean action |jTo|, |ll| 



i-np rl • e 1 
5[R(r)] = / dr -mB} + i-{B^x - B\\z)y + V{z) + -mul^x^ + y"^ 
JO L2 c 2 
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where (3 = l/ksT is the inverse temperature. The tunneling rate is proportional to exp{—Sc/fi}, 
where the semiclassical action Sc is determined by the bounce solution of the equation 6S = 0, 
in which the periodic boundary condition, R(^/3) = R(0), is required. This is so-called the 
method of the imaginary part of the free energy F = — fc^TlnZ |T2| calculated from the 
partition function 

Z = J dK' D(R';R') , (4) 

with the density matrix 

D{R'{0);R"{hp)) = Jl^'R exp|-i5[R]| . (5) 

This method is identical to the WKB method at zero temperature, and allows us to have a 
unified treatment of the escape rate for finite temperatures. 

We are interested in the electron tunneling out of the metastable state z = in the 
z-direction. After the tunneling other degrees of freedoms, x and y, can take any allowed 
value. Therefore the summation over final states, integrations over the x and y coordinates, 
will be taken in the calculation of the partition function. Those are gaussian integrals. To 
integrate over the x-coordinate, we perform a Fourier transformation on the time interval 

[0, hpy. 

-| oo 

(a;(r),2/(r),z(r)) = — ^ (x., • (6) 
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Here = 2Tm/hf3 . The action 5* of eq.(3) can be then rewritten as 



1 oo 
n=-o< 



]^mvlznZ-n + V{z)n + ]j{invl + mujl)xnX. 



1 °° rl . e 

+ ^ ni^^l + mU}o)yny-n + -{B^X^ - 5|| 

'^P n=-oo 



(7) 



Here Un — —y-n has been used. Now, we can shift the origin of {yn}, and integrate over 
them. Then we obtain the effective action as 



'^P n=-c 



1 ^ n 1 



1 ~ 1 



Similarly, the integration over can be performed. The resulting effective action is 



[z{t)] = / dr 
Jo 



-m^ + V{z) 
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with the first kernel 



and the second kernel 



- 92{r - t')][z{t) - z{t' 



1 °° 1 i^BwYv^ 



(9) 



(10) 
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to 
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n=-oo ^ [777,^/2 _^ ma;o] ( ) + ("^^n + ?^'^o)^ 

Here /J*^ dr ^i(r) = and dr 9^2 (t) — have been used in arriving above result. 

Eqs.(9-ll) can be simplified. We note that the kernel can be separated into two parts: 



with 
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The first part of the right hand of eq.(12) is simply a periodic delta function, which gives no 
contribution to the effective action and will be dropped. We further simplify the expressions 
by taking the limit of cuq — 0, and find that 



9i{r) = -9i{t) = - 



1 m feB\\ 



(14) 



and 



m f eB\\ 
2 V mc 



sinh 



hp 



(15) 



with LJ± = B^e/mc . Then using eqs. (14,15,9) we obtain the effective action for our problem 



as 



hf3 



dr 



-mz'^ + V{z) 



2 Jo Jo 2 \ mc 



00 ± cosh[a;j_( 
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r — r 



sinhf^ 



[z{t) - Z{' 



'M2 



(16) 



Thus, we have mapped the original 3-d problem onto a 1-d one. This is because the in- 



plane motion, x- and ?/-directions, effectively behaves as a bath for 2;-direction motion [pi) 
where the parallel magnetic field B\\ serves as a coupling coefficient. This form of effective 
action is easy to handle because of the existence of a sophisticated technique to calculate the 
dissipative tunneling rate|TO[. In terms of the dissipative tunneling[p!0[], the environmental 
effect in the eq.(16) is a superohmic damping of s = oo for Bj_ ^ 0, where the environment 
consists of the in-plane cyclotron motion. If B± = 0, it is then a subohmic damping of 
s = 0, where the environment consists of in-plane plane-waves. We should point out that 
no approximation is used in obtaining eq.(16). Simply by inspection of eq.(16), one can 
conclude that the effective bath due to the in-plane motion will influence the tunneling rate 
out of the confining plane by affecting the semiclassical action. A numerical calculation of 
the tunneling rate will be needed to cover whole parameter region. In the following we take 
B_i = 0, and focus on the peculiar subohmic damping case with B\\ ^ 0. 
From eq.(16), the effective action with B± = is 



SefMr)] = dr 
^ u 



1 .2 

-mz 
2 



Viz) 



m I eB\\ 



hfi 



1 r^P 

. dr dT'[z{r) - ziT')f . (17) 
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This equation is explicitly gauge invariant under the change y ^ y + constant, because of 
the periodic boundary condition of x, y imposed in the tunneling calculation as pointed out 
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above. If there is a tunneling solution, it can be shown that 



lim / dr z{t) 



constant 



for the semiclassical solution. This is the result of the period boundary condition that 
lim;j^_^oo z{hf3) = z{0) = and the finiteness of the semiclassical action. Therefore, the 
effective action at zero temperature may be written as 



Seff[z{r)] 



dr 



-mz\T) + V{z{T)) + ji — 



zHr) 



mc 



(19) 



It shows that the original potential V{z) is renormalized to V{z) + m (ei?||/mc)^z^/2 , 
and the metastable state at 2; = becomes more stable. For a large z, the renormalized 
potential is positive now, and we have a double- well like potential |jl3|. If the second local 
potential minimum is lower than the one at 2; = 0, tunneling out of 2; = is finite. However, 
for a large enough parallel magnetic field -By, the second local potential minumum will be 
higher than the one at z = 0. The tunneling rate is then zero. In this case, eq.(18) is not 
valid. This analysis suggests that there is a magnetic field induced localization transition at 
zero temperature. The violation of eq.(18) may serve as the indication for this transition, or 
alternatively, the comparison of the local two potential minima of the renormalized potential 
in eq.(19) may determine the critical magnetic field. Thus we have obtained that the s = 
dissipative environment is marginal for localization in tunneling from a metastable state. 



compared to the s = 1 case for the tunneling splitting||14|, |15 



Away from the localization region, for a small parallel magnetic field Bi^\ the tunneling 
rate is finite. In particular, for a very small parallel magnetic field and low temperatures the 
semiclassical action may be evaluated perturbatively: 



dr 



mziiT) + Viz,iT)) 
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m 



'eBii 



knT 
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dTzJr] 
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where Zc{t) is the bounce solution at zero temperature without the magnetic field, determined 
by the equation = V{zc). Eq.(20) shows a pronounced i?y dependence and linear 

temperature dependence. 

Now we briefiy discuss the effect of finite perpendicular magnetic field Bj^. For ^ 



the effective bath changes from s = at = to s = oo. From Refs. |[T5|, we expect 
that the tunneling can occur at any value of B\\ and there is no localization. Therefore by 
tilting the direction of the magnetic field there should be a large change in the tunneling 
rate. 

To conclude, we discsuss the experimental verification of the above results. The tunneling 
rate of trapped electrons escaping from a helium- vacuum interface was measured in Ref . @] as 
a function of -Bjj and T with B± = 0. A numerical estimation shows the data are in the low 
parallel magnetic field and low temperature limit. The comparison between the calculation 
according to eq.(20) and the data of Ref.[^ is shown in the Figure. A quantitative agreement 
is found in the low temperature limit. Note that there is no fitting parameter. The deviation 
in the high temperature end in the Figure suggests that the thermal activation starts to play 
a role. It is interesting to note that the corresponding critical magnetic field according to 
the renormalized potential in eq.(19) is about 1.3 Tesla for that experiment, which seems 
easy to realize experimentally. In Ref.|Q it is also found that there is no dependence of 
tunneling rate on 5^ when B\\ = 0. This is implied in eq.(16), too, where by taking B\\ = 
the in-plane motion is decoupled from the tunneling in z-direction. To further assure oneself 
of the agreement between eq.(20) and the experimental data, the absence of the effect of the 
other dynamical correlations in the experiment need to be addressed. We discuss this below. 

It has been shown that the response of the remaining 2-d electrons to the tunneling elec- 
tron is superohmic0. The infiuence from ripplons is not only superohmic, also weak. The 
only possible ohmic effect is from the He atom scatterings. However, its strength is expo- 



nentially small at low temperature, and can be ignored at low temperatures. Because the 
superohmic bath has a weak influence on tunneling in the absence of a magnetic field[|l^, |12|, 
the tunneling is therefore dominated by the adiabatic potential. Therefore the adiabatic po- 
tential is good in the explanation of the experiments |1|] and dynamical correlations from 
those environments are weakQ. Furthermore, it has been shown recently that, the super- 
ohmic bath has a weak influence on the tunneling rate even in the presence of a magnetic 
field [0. Therefore the adiabatic potential can be used in the tunneling rate calculation with 
a magnetic field, too. Using the adiabatic potential first proposed by lye et al in Ref.|^, we 
have calculated the semiclassical action by eq.(20) as shown in the Figure. In the calculation 
the Stark shift and the distance of the 2-d electron layer from the helium-vacuum interface 
have been considered. 
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Figure Caption. 

Figure. The circles are the experimental data from Ref.[3]. The solid line is the calculation 
according to eq.(20), in which the 2-d electron density n = 8.7 x and the external 

electric field E — 30V/ cm are used. 
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